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Abstract 

In my lectures I shall give an introduction to the "membrane model" of 
black holes (BHs). This is not a new theory, but a useful reformulation of the 
standard relativistic theory of BH -as far as physics outside the horizon is 
concerned- which is much closer to the intuition we have gained from other 
fields of physics. While the basic equations look less elegant, it has the advantage 
that we can understand astrophysical processes near a BH much more easily. 
In the membrane model one first splits the elegant 4-dimensional physical laws 
into space and time (3+1 splitting). For a stationary BH there is a preferred 
decomposition. Relative to this the dynamical variables (electromagnetic fields, 
etc.) become quantities of an "absolute space" which evolve as functions of 
an "absolute time", as we are accustomed to from nonrelativistics physics. We 
shall see, for example, that the 3+1 splitting brings Maxwell's equations into 
a form which resembles the familiar form of Maxwell's equations for moving 
conductors. We can then use the pictures and the experience from ordinary 
electrodynamics . 

In a second step, one replaces the boundary conditions at the horizon by 
physical properties (electric conductivity, etc.) of a fictitious membrane. This 
procedure is completely adequate as long as one is not interested in fine de- 
tails very close to the horizon. The details of this boundary layer are, however, 
completely irrelevant for astrophysical applications. 

The following points will be discussed in detail: 

• Solutions of Maxwell's equations in a Kerr background 

• Space-time splittings 

• The horizon as a conducting membrane 

• Magnetic energy extraction from a BH 

• BH's as current generators or rotators of electric motors 

• The Blandford-Znajek process 

• Stationary axisymmetric electrodynamics for force-free fields 

1 Introduction 

In these lectures I give an introduction to what is called the membrane model of black 
holes (BHs). This is not a new theory, but a convenient reformulation of the standard 
relativistic theory of BHs - as far as physics outside the horizon is concerned -, which 



is much closer to the intuition we have gained from other fields of physics (see Q, for 
a general reference). While the basic equations look less elegant, it has the advantage 
that we can understand astrophysical processes near a BH much more easily. For 
an analogy, imagine you would have to explain how a Tokomak works by using the 
language and pictures of special relativity (SR), i.e., by using the electromagnetic field 
tensor and 4-dimensional pictures of plasma flows. I would not know how to do this and 
how to get, for instance, an understanding of even the simplest plasma instabilities. A 
closer analogy would be to translate relevant studies of the electrodynamics of pulsars 
into a 4-dimensional language. The basic equations look beautiful, but it would be hard 
to understand anything. (You may say that radio pulsars are anyhow not understood.) 

In the membrane model (often called membrane paradigm 0) one first splits the 
elegant 4-dimensional physical laws of general relativity (GR) into space and time 
(3+1 splitting). For a general situation this can be done in many ways (reflecting 
the gauge freedom in GR) since there is no canonical fibration of spacetime by level 
surfaces of constant time. However, for a stationary BH there is a preferred decom- 
position. Relative to this the dynamical variables (electromagnetic fields, etc) become 
quantities on an absolute space which evolve as functions of an absolute time, as we 
are accustomed to from nonrelativistic physics. We shall see, for example, that the 
3+1 splitting brings Maxwell's equations into a form which resembles the familiar 
form of Maxwell's equations for moving conductors. We can then use the pictures and 
the experience from ordinary electrodynamics. 

In a second step one replaces the boundary conditions at the horizon by physi- 
cal properties (electric conductivity, etc) of a fictitious membrane. This procedure is 
completely adequate as long as one is not interested in fine details very close to the 
horizon. The details of this boundary layer are, however, completely irrelevant for 
astrophysical applications. (The situation is similar to many problems in electrody- 
namics, where one replaces the real surface properties of a conductor and other media 
by idealized boundary conditions.) 

The program of these lectures is as follows. First I will discuss the 3+1 splitting 
of the spacetime of a stationary rotating BH and of Maxwell's equations outside its 
horizon. We shall see that this can be achieved very smoothly by using the calculus 
of differential forms. As an illustration and for later use we shall apply these tools 
for a discussion of an exact solution of Maxwell's equation on a Kerr background, 
which describes an asymptotically homogeneous magnetic field. We shall then derive 
the electromagnetic properties of the fictitious membrane that simulate the boundary 
conditions at the horizon. Here, I can offer a much simpler derivation than has been 
given so far in the literature. As an important example of a physical process rela- 
tively close to a BH I will treat in detail the magnetic energy extraction of a hole's 
rotational energy. Blandford and Znajek have first pointed out the possible relevance 
of this mechanism for an understanding of active galactic nuclei. It may well play 
an important role in the formation of energetic jets. The Blandford-Znajek process 
could also be important for explaining gamma-ray-bursts, because it may energize a 
Poynting-dominated outflow. 

I hope to show you that the physics involved is not very different from that behind 
the electric generator in Fig. [I]. 




Figure V. Electric generator whose physics is similar to the electrodynamics of black 
holes in external magnetic fields. 

2 Space-Time Splitting of Electrodynamics 

I describe now the 3+1 splitting of the general relativistic Maxwell equations on a 
stationary spacetime (M,^g). Most of what follows could easily be generalized to 
spacetimes which admit a foliation by spacelike hypersurfaces (see, e.g., Ref. 0), but 
this is not needed in what follows. 

Slightly more specifically, we shall assume that globally M is a product R x E, 
such that the natural coordinate t of R is adapted to the Killing field k, i.e., k = dt- 
We decompose the Killing field into normal and parallel components relative to the 
"absolute space" (E, g), g being the induced metric on E, 

d t = au + /3. (1) 

Here u is the unit normal field and (3 is tangent to E. This is what one calls the 
decomposition into lapse and shift; a is the lapse function und (3 the shift vector 
field. We shall usually work with adapted coordinates = (t,x l ), where {x 1 } is 
a coordinate system on E. Let (3 = (3 l di {pi = d/dx l ), and consider the basis of 
1-forms 

a dt, dx l + (3 { dt. (2) 

One verifies immediately, that this is dual to the basis {u, <9,} of vector fields. Since u 
is perpendicular to the tangent vectors <9j of E, the 4-metric has the form 



W g = - a 2 dt 2 + g ti (dx i + ftdt) (dx j + (3 j dt) , 



(3) 



where gijdx l dx^ is the induced metric g on E. Clearly, a, /3, and g are all time- 
independent quantities on S. 

For what follows, I would like to change this setup slightly by using, instead of di 
and dx l , a dual orthonormal pair {e^} and on E. Instead of @, we have then the 
orthonormal tetrad 

0° = adt, 0* = 0* + l?dt, (4) 
where now (3 = f3 l ei. This is dual to the orthonormal frame 

e = u = —(dt - 0), ej. (5) 
a 

The tetrad {0^} describes the reference frames of so-called FIDOs, for fiducial ob- 
servers. Their 4- velocity is thus perpendicular to the absolute space E. 

Relative to these observers we have for the electromagnetic field tensor F (2-form) 
the same decomposition as in SR: 

F = E A 6° + B, (6) 

where E is the electric 1-form E = E;fi l and B the mag netic 2-form B = \B^ 6 l A Q 3 . 
(Ei, Bij are the field strengths measured by the FIDOs.) 

In a second step we decompose E and B relative to absolute space and absolute 
time. We have, using ([5]), 



where 



Similarly, 



E = E i 6 i = Ei (0* + pdt) =£ + tp£dt, (7) 



£ = Erf 1 , i, : interior product. (8) 



B = B + dt AipB, £ = # A (9) 



Together we arrive at the following 3+1 decomposition of F: 

F = B+ (a£ -ipB) Adt. (10) 

^From this the 3+1 splitting of the homogeneous Maxwell equations is readily 
obtained: dF = gives 

dB + dt A d t B + d(a £) Adt- B) A dt = 0. 

Here d denotes the differential on E. This gives the two equations 

dB = 0, d(a£) + d t B = d(ipB) 

or, with the Cartan identity Lp = d o ip + o d, 

di3 = 0, d(a£) + (d t — Lp)B = 0. (11) 



The second equation describes Faraday's induction law in a gravitational field. It will 
be of crucial importance in later sections. Note, in particular, the coupling of the 
£>-field to the shift through the Lie derivative. 

Let us also decompose the representation of F by a potential, F = dA. We have, 
using again (^J), 

A = A^e^ = aA G dt + A i (d i + (3 i dt) 
= (aA + i(3 A)dt + A, A = A i 'd i . 

Thus 

A=-(j)dt + A, (12) 



where 



This gives 



■(aAo + ifsA). (13) 



dA = -d(j) Adt + dA + dtAdt A, 

which is of the form (|10D , with 

B = dA, aS = -d<p- dtA + ipdA. (14) 

Apart from the last term, this is what one is used to. 
Now, we turn to the inhomogeneous Maxwell equation 

d*F = 4nS. (15) 

We need first the Hodge-dual of (JTTJ|) . We decompose *F similarly to ([]): 

*F = -H A 6° + D, (16) 

which can be viewed as a definition of H and D. Comparison with (|j) shows, that 

H = BiP, B 1 = B 23 , etc, 

D = Ei 8 2 A6 3 + E 2 6 3 A8 1 + E 3 8 1 A 9 2 . (17) 
With (^) we find (* denotes the Hodge-dual on E) 

H = H + ipHAdt, H = *B, 

D = V-ipVAdt, V=*S. (18) 
If this is inserted into fliTf), we obtain 

*F = V - (aH + ip V) A dt. (19) 
The dual J = *S of the current 3-form can be decomposed as in SR 

J = Pel0°+Jk0 k , (20) 



where p e i is the electric charge density and j k is the electric current density relative to 
the FIDOs. We use them to introduce the following quantities on the absolute space 

p = p el $ 1 A$ 2 A$\ j = j k # k , J = *3- (21) 
Using the notation 77^ := *6 fl , we can decompose S as follows 

S = pelV°+JkV k 

= Pel e 1 a e 2 a e 3 - & e 2 a e 3 + . . . ) a e° 

= p + peiifi 1 tf 2 A t? 3 + . . . ) A dt - a (ji ^ 2 A i? 3 + . . . ) A dt. 

Thus 

S = p+(ipp-aj) Adt. (22) 

Inserting this and (|19|) into ([15]) leads to 

d*F = dV + dtAd t V-d(aH)Adt-d(ip'D)Adt 
= Anp + 4n(ip p - a J) A dt, 

and hence to the following 3+1 split of the inhomogeneous Maxwell equation 

dV = Anp, d{aH) = {d t - Lp)V + 4vra J. (23) 

^From these laws one obtains immediately the local conservation law of the electric 
charge (use that d commutes with Lp): 

(d t -Lp)p + d(aJ) = 0. (24) 

This follows, of course, also from dS = and the decomposition (|22j). 

Integral Formulas 

As is well-known from ordinary electrodynamics, it is often useful to write the 
basic laws (TT|), (p3|), and (|24| ) in integral forms. Consider, for instance, the induction 



law in (|Tl ) . If we integrate this over a surface area A, which is at rest relative to the 



absolute space, we obtain with Stokes' theorem (C := OA) 

a£ = -^ + ! B+ ! Lp B. 
c dt J A J A 

Here, we use LpB = dipB (since dB = 0) and Stokes' theorem once more, with the 
result 



f d 
at = — — 

dt 



B + jipB. (25) 



The left hand side is the electromotive force (EMF) along C. The last term is similar to 
the additional term one encounters in Faraday's induction law for moving conductors. 
It is an expression of the coupling of B to the gravitomagnetic field and plays a crucial 



role in much that follows. This term contributes also for a stationary situation, for 
which (EH) reduces to 



EMF(C) = j aS = j ipB. (26) 

The integral form of the Ampere-Maxwell law is obtained similarly. Integrating 
the second equation in (EjJ), we obtain with the Cartan identity L^do^ + i^od 
and Gauss' law (first equation in (p3|)): 

j{aH + ipV) = -J^D + A-K J {aj-ipp). (27) 



The integral form of charge conservation is obtained by integrating (p4|) over a 
volume V which is at rest relative to absolute space: 

d_ 
dt 



p= - / {aj-ipp) (28) 

V JdV 



(note that Lp p = dip p). 

One could, of course, also derive integral formulas for moving volumes and surface 
areas (exercise). 

Vector Analytic Formulation 

The similarity of the basic laws in the 3+1 split with ordinary electrodynamics 
becomes even closer if we write everything in vector analytic form. I give a dictionary 
between the two formulations that is valid for any 3-dimensional Riemannian manifold 

(£,*)■ 

The metric g defines natural isomorphisms Jj and b between the sets of 1-forms, 
A X (E), and vector fields, X(Y0). In addition, the volume form 77, belonging to the 
metric, defines an isomorphism between X(H) and the space of 2-forms, A 2 (E), given 
by 

B\ — >B = i s r). (29) 

We have the following commutative diagram, in which * denotes, as always, the Hodge- 
dual: 

A X (S) * ! A 2 (E) 





^From this one can read off, for instance, 

izr} = *v (v e X(Y), v : (vf). (30) 
The cross product and the wedge product are related as follows: 

A X (S) x A : (S) a 2 (E) 

1 1 
X(E) x X(E) X(E) 



In particular, we have 

ivxwV = vAw - (31) 

With the help of the next commutative diagram one can reduce many of the vector 
analytic identities to d a d — 0. 

— ► A°(E) -±> A X (E) -±> A 2 (E) A 3 (E) — > 

II b tl tl 1 i.v T •»? 

_> A°(S) ^ #(£) ^ AT(E) ^> A°(£) — > 0. 

We can read off, for example, 

icmiffV = dv- (32) 
For a 1-form w with vector field w = w^,we have the translation 



i$dw icuxwicmiwV = [(curlw) x w] b . (33) 
Here, we made use of the algebraic relation 

ifficrj ® z* * u = *(u A v) ® (u x tT) b . (34) 
We need also the translation of the Lie derivative of a 1-form u>: 



Ltfiy = d itfw +irfdw — d(v,w) + [(curlw) x vf. 



(v,w) 

Thus, 

L$w = {grad (u, w) + (curliu) x v} h . (35) 
Similarly, we have for a 2-form B = isT]\ 

L^B = L$ ig rj = i$ d ig rj + d i$ ig rj 

divBr) (Bxv)^ 

@ . 

~~ l {(divB)v + curl(Bxv)} r l- 

We have thus the correspondence 

L#B < — > (divB)v + curl (B x v). (36) 

Here we have to stress that the right hand side is in general not equal to L$ B = [v, B] 
(Lie bracket). This comes out as follows: 

L$ B = Ltfigri = [La, ig] rj + ig Lffj 

= l {(diYv)B + [v,B]} T 1- 

The correspondence ([36|) is thus equivalent to 

Li? B< — ► Li?B+ (div v) B. (37) 



— * 

Only for divv = do the Lie derivatives L^B and L^B correspond to each other! 
In Maxwell's equations the Lie derivative L^B occurs. This will be replaced in 

the vector analytic translation by L^B, because div/3 is (for a stationary metric) 
proportional to the trace of the second fundamental form of the time slices and this 
vanishes for maximal slicing. For the Kerr solution we are, for instance, in this situation 
(exercise). 

Part of what has been said is summarized for convenience in the table below. 

Dictionary 



calculus of forms 


vector analysis 


notation 


v Aw 


V x w 


v = (vY, w = (w) b 




B x v, E x v 


B = i S Tj, V = igrj 


df 


grad/ 


f : function 


dv 


curl {7 




dB, dV 


divB, div E 






grad (v, w) — v x curl w 


w = (wY 


Ly B , Ly D 


(div B)v — curl(-i7 x B), B < — > E 





Summary 



For reference, we write down once more the 3+1 split of Maxwell's equations fllTD 
and (^3|) in Cartan's calculus 

dB = 0, d(a E) + (d t - L P )B = 0, 
dV = 4vrp, d(aH) = {d t - Lp)V + Ana J. (38) 

The dictionary above allows us to translate these into the vector analytic form: 

V-S = 0, V x (aE) + (d t -L$)B = 0, 
V-E = 4irp el , V x (aB) = (d t - L^)E + 4vraj. (39) 

3 Black Hole in a Homogeneous Magnetic Field 

As an instructive example and a useful tool we discuss now an exact solution of 
Maxwell's equations in the Kerr metric, which becomes asymptotically a homogeneous 
magnetic field. This solution can be found in a strikingly simple manner ||. 
For any Killing field K one has the following identity 

5dK b = 2R(K), (40) 

where 5 denotes the codifferential and R(K) is the 1-form with components R 11V K V . 
In components ( [40|) is equivalent to 

K» a , a = -R» a K a . (41) 

This form can be obtained by contracting the indices a and p in the following general 
equation for a vector field 



and by using the consequence K a a = of the Killing equation K a . p + K p . a = 0. 
For a vacuum spacetime we thus have 

8dK b = (42) 

for any Killing field. Hence, the vacuum Maxwell equations are satisfied if F is a 
constant linear combination of the differential of Killing fields (their duals, to be 
precise). For the Kerr metric, as for any axially symmetric stationary spacetime, we 
have two Killing fields k and m, say; in adapted coordinates these are k = dt and 
m = dp. The Komar formulae provide convenient expressions for the total mass M 
and the total angular momentum J of the Kerr BH: 

M = -—[ *dk\ J=—[ *dm b (43) 
8tt 16tt 

(for G — 1). 

We try the ansatz 

F = ~ B (dm? + 2a dk b ) (B = const), (44) 
and choose a such that the total electric charge 



47T 

vanishes. The Komar formulae (E3f) tell us that 



Q = — — — / *F (45) 



Q = — - 5 Q (16tt J - 2a • 8ttM), (46) 
8n 

and this vanishes if a = J/M (which is the standard meaning of the symbol a in the 
Kerr solution). 

Clearly, F is stationary and axisymmetric: 

L k F = L m F = 0, (47) 

because (dropping b from now on) 

Lf~dk = d Lk k = (L^ A; = [k, k] =0), etc. 

For the further discussion we need the Kerr metric. In Boyer-Lindquist coordinates 
and more or less standard notation it has the form (^), i.e., 

{4) g = [-a 2 dt 2 + g w (d V + f3 v dt) 2 } + [g rr dr 2 + g M d$% (48) 

with only the component /3 V of the shift being ^ 0. With the abbreviations 

p 2 ■= r 2 + a 2 cos¥, A := r 2 - 2Mr + a 2 , S 2 := (r 2 + a 2 ) 2 - a 2 A sin 2 $, (49) 

the metric coefficients are 

n 2 Y 2 

P 2 • 2 q 

9rr = x, 9m = P , SW = sm ^ -T> 
A 

2Mr 2Mrasin 2 i? 

= -H 5~, ^ = 5 ' ( 50 ) 

p z 



while the lapse and shift are given by 

2 P 2 A 

This gives asymptotically 
2M 



2Mr 



(51) 



(i) 



9 



dt 2 



4aM 



sin't? + O ( — 



dt dip 



+ 



l + O ( - 

r 



[dr 2 + r 2 (d$ 2 + sin 2 $d<^ 



(52) 



To establish the connection with our general discussion, we introduce here as orthonor- 
mal basis of the absolute space naturally 



$ r = sfg^-dr, ^ = y/gmdfi, = y/g^dip. 

The shift vector is (3 = fi* d v . 

The angular velocity oj of the FIDOs, with 4- velocity u = -{d t — /3 V d v ), is 



(53) 



u 



W = -r = -P 



ir 



9up_ 



(54) 



The last equality sign implies that the angular momentum of these (Bardeen) observers 
vanishes: (u, m) = 0. The FIDOs are, therefore, sometimes also called ZAMOs (for 
zero angular momentum observers) . 

Now, we want to discuss in detail the solution fl44]), which we can express in terms 
of a potential: F = dA, with 



A = - B (m + 2ak). 



(55) 



Let us first look at the asymptotics. The 1-forms k^dx 11 , m fJi dx fl belonging to the 
Killing fields (k^ = g^t, m M = g^) are asymptotically k ~ — dt, m ~ r 2 sin 2, # dip, 
whence ( fl"4"D gives 

F ~ Bq [sin $ dr A r sin dip + cos d rdd A r sin d dip] . (56) 

This is a magnetic field in the z-direction whose magnitude is Bq. 
In fl5"5|) we need 

m + 2ak = (p w + 2a o Mt ) dx M = (g^, + 2a gr tt ) dt + (g w + 2a 5^) 



'-ojg w + 2ag tt )dt + (g w - 2auj g w )dip. 



Using the notation 



~ 2 _ 



u) = — sin-i? 



we thus have 



m 



+ 2a£; = [-wcj 2 + 2a (, 



2-2 
UJ UJ 



a 2 )] dt + uj 2 {l -2auj)dip. 



(57) 



Comparing this with ((T^j), we obtain for the potentials 

1 



A 



B [uu 2 + 2a (a 2 -uj 2 uj 2 )] 



A^dif, A<p = - B ti 2 (l - 2au). 



The fields £ and B can now be obtained from (0). We have 
B = dA = Aip >r dr A dip + A v ^ d$ A d<p 

A v o f A & 3 -VAA<p r f A -&\ 



1 



E sin-$ 



.4.^ is explicitly (use fl59l), (|57|) , (154]), and (|51 



^= - S — sin 2 ?? ( I - if/ 



Mr 



We write this as 



A v — -BqX, 



p 



With this notation, ( |5T!D reads 



B 



Bn 



2Esin$ 

The corresponding vector field B is thus 

B 



B 



X # e r — v A X r e$ 



2Esin$ 

For £ we have, with (3 = —uj d^, 

a £ = —dcj) + ip dA = —d(f> + uj dA v . 
[From this one finds quickly 



E 



+ 



S aE f 


"9a 2 


P 2 


1 


dr 


1 


r&* 2 

h r 



Msin 2 $ 
^ 2_ 



E 2 - 4a 2 Mr) 



_9_ 



(58) 
(59) 



(60) 



(61) 

(62) 

(63) 
(64) 



(65) 



The field lines of E are shown in Fig. 

It is of interest to work out the magnetic flux through the equator of the BH, i.e., 



B 



upper h. 



A = 2tiA, 



ip | equator 



(66) 



equator 



Specializing (|6"ID to d = n/2 and r = rn gives (the horizon is located at A = 0) 

$ = AnBo M(r H — M) — AnB M^M 2 - a 2 . (67) 

Here, we have used 



r H = M + ^fW 



a". 



(6* 



Note that this vanishes for an extremal BH (a = M). Generically one has, as expected, 



Figure 2: Field lines of E (from Ref. [1]). 



4 The Horizon as a Conducting Membrane 

As long as one is not interested in fine details very close to the horizon, one can regard 
the boundary conditions implied by the horizon as arising from physical properties of 
a fictitious membrane. This membrane is thought of as endowed with surface charge 
density cr#, surface current Jh-, and surface resistivity Rh- This idea was first pursued 
by Damour M and independently by Znajek 0. Later it was further developed by 
Thorne and Mac Donald || and other authors (see |J and references therein). 

Below we give a much simplified derivation that the following surface properties 
for the parallel and perpendicular components hold: 



Gauss's law 
Ampere's law 
charge conservation 
Ohm's law 



E± — >Akg h , 

a B\\ — > B H = 47r Ju x n, 

aj± — > -d t a H - (2) V -J H , 

a En — >E H = R H J H . (69) 



The surface resistivity Rh turns out to be equal to the vacuum impedance 

R H = 4tt = 377 Ohm. (70) 

For the horizon fields Eh, Rh we have, therefore, 

Bh = E H x n, (71) 

as for a plane wave in vacuum. 

Toward the horizon the FIDOs move relative to freely falling observers, say, more 
and more rapidly, approaching the velocity of light. Mathematically, the tetrad {0^} 



becomes singular at the horizon, and therefore the components of F relative to {9^} 
are ill-behaved. The 2-form F should, of course, remain regular and the laws ( |SUD 
and fl70|) are just an expression of this requirement. For illustration, we demonstrate 
this first for a Schwarzschild BH, and generalize afterwards the argument in a simple 
manner. 



(a) Derivation for a Schwarzschild BH 

The procedure is simple: we pass to a regular, not necessarily orthonormal tetrad. 
The angular part 9 2 = rdf}, 9 3 = rsmfidtp is kept, but instead of 9° = adt, 9 1 = 
oT x dr we use di, where t is the Eddington-Finkelstein time coordinate 

t = t + 2M\n[—-l) — >dt = dt + a~ 2 — dr (r > 2M). (72) 

Since the Eddington-Finkelstein coordinates are regular at the horizon, the same is 
true for the basis {di, dr, 9 2 , 9 3 }. We have the relations 

0° = adt dr, 

a r 

9 1 = -dr, (73) 

a 

which allow us to rewrite the decomposition (^j) as follows 

F = E 1 9 1 A9° + + B^ 1 A9 2 + ... 

= E 1 dr Adi+ E 2 (a9 2 Adt-- - — 9 2 A dr) + E 3 
\ a r 

+B 3 - dr A 9 2 - B 2 - dr A 9 3 + B 1 9 2 A 9 3 

a a 



or 



E x dr A dt + a E 2 9 2 A dt + a E 3 9 3 A dt + B x 9 2 A 9' 3 
+±(b 3 + —E 2 ) drA9 2 + -(-B 2 + —E 3 ) drA 9 3 . (74) 



a \ r j a \ r 

The regularity of the coefficients in this expansion implies the following behavior when 
the horizon is approached (a | 0): 

(z) radial components: E\, B\ = 0(1), 
(ii) tangent components: E 2 , E 3 , B 2 , B 3 = O ( — 
(ill) E 2 + B 3 ,E 3 -B 2 =0(a). (75) 



This shows that 



aE\\ = n x aB\\ + 0(a ), (n:=e r ), 
aB\\ = -fix aE\\ + 0(a 2 ), (76) 
E n = E±, B n = Bj_ remain finite. (77) 



This is basically already what was claimed in ( |69j) and ( ffD|) . It is useful to introduce 
the stretched horizon 7i s = {a = an C 1} which is arbitrary close to the event 
horizon. Let 



^ : =(«£||) QH , B H :={aB\\) aH 

and as above 

— * — * 

E n := E ■ n, B n := B ■ n. 
These components remain finite for a fl { and we have up to 0(a H ) 

E H = n x B H , Bh = —n x E H . 



(7i 



(79) 



JO) 



The surface charge density an and the surface current density Jh on 7i s are defined 
by 



En\ 

4?r y 



B 



H 



An J H xn 



*1) 



47T. 



The second equation and (|q!]| ) imply Ohm's law in 

Finally, we make use of the Ampere-Maxwell law (|39|) (for (3 = 0): 

d t E = V x (afl) - 4vra J. 
The normal component on T~C S is 



*2) 



V x (aB) 



An a j n 



H 3 



Using 



Vx(tt5) 



-4tt {2) V-J h , 



V x (4tt Jj/ x ft) 

where denotes the induced covariant derivation on H, s , we obtain 

d t o- H + {2) V ■ J H + (aj n ) n s = 0, 



13) 



as an expression of charge conservation. 

This completes the derivation of ( |69"D and (|70|) for the Schwarzschild BH. Next, we 
generalize the discussion to an arbitrary static BH. 

(b) Derivation for static BH 

We introduce first a parametrization of the exterior metric which was used also in 
Israel's famous proof of the uniqueness theorem for the Schwarzschild BH. 
The starting point is (0) for (3 = 0, i.e., 



( 4 ), 



g = -a 2 dt 2 + g, a and g independent of t. 



*4) 



Note that a 2 = —(k\k) > 0, k = d t , and that the horizon has to be at a = 0. We 
assume that the lapse function has no critical point, da 7^ 0. The absolute space E is 
then foliated by the leaves {a = const}. The function 

p := (da\da)~^ (85) 

is then positive on E. 

Now we introduce adapted coordinates on E. Consider in any point p G E the 
1-dimensional subspace of T P (E) perpendicular to the tangent space of the leave 
{a = const} through p. This defines a 1-dimensional distribution which is, of course, 
involutive (integrable). The Frobenius theorem then tells us that we can introduce 
coordinates {x 1 } on E, such that x A (A = 2, 3) are constant along the integral curves 
of the distribution. For x 1 we can choose the lapse function, and thus obtain 

g = p 2 da 2 +g, g = cjABdx A dx B . (86) 

Here, p and qab depend in general on all three coordinates x 1 = a, x A . 

We also need the surface gravity k on the horizon. A useful formula is (see, e.g., 

0) 

K 2 = -hdk\dk)\ R . (87) 

Now, k = —a 2 dt, dk = —2a da A dt, whence 

k = — . (88) 
Ph 

^From the zeroth law of BH physics we know that k = const. Since we want to assume 
a regular Killing horizon (generated by k), we conclude 

< Ph < oo, Ph = const. (89) 



Combining (84) with (jSq) we have outside the horizon 



2 



^g = -Ndt 2 + ^dN 2 +g, N:=a 2 . (90) 



The natural FIDO tetrad is 



6° = VNdt, 

9 1 = ^=dN, 
2VN 



9 A (A = 2, 3) : orthonormal 2-bein for g. (91) 

This becomes again singular at the horizon (N = 0). 

Now we imitate what we did for the Schwarzschild BH. We search for a basis of 
1-forms which is well-defined in the neighborhood of the horizon. Guided by fl?2| ) we 
introduce 



9 t = dt + — dN 
2N 



(92) 



and rewrite (|9CiD 



(4) 



g = -N(e t ) 2 + p9 t dN + g. 



Thanks to (p9|) we conclude that 

{9\ dN,9 A (A = 2,3)} (93) 

remains a regular basis on the horizon. Outside the horizon we can express (|91|) in 
terms of this basis: 

6° = VN6 1 ?=dN, 

9 1 = -?=dN. (94) 
2^N 



We can now proceed as in the derivation of ([71]), obtaining now 
o E 

F = f-^dNA6 t + VNE 2 9 2 A9 t + VNE 3 e 3 A9 t + B 1 9 2 A9 3 

+ (E 2 + B 3 ) dNA9 2 + -B— (B 2 - E 3 ) dN A 9 3 . (95) 
2yN 2\/N 



This implies again the limiting behavior ( fT5|) for the normal and parallel components 
of the electric and magnetic fields. 

(c) Derivation for rotating BHs 

Finally, I give a similar derivation for rotating BHs. (This was worked out in 
collaboration with Gerold Betschart in the course of his diploma work.) 

I first need the Papapetrou parametrization of a stationary, axisymmetric BH. 
Since this will be treated in detail in the lectures by Markus Heusler 0, I can be 
brief. 

Since the isometry group is R x 5*0(2) we have two commuting Killing fields k 
and m, say, which are tangent to the orbits belonging to the group action. We assume 
that k and m satisfy the Frobenius integrability conditions 

k A m A dk = 0, k A m A dm = 0. (96) 

The Frobenius theorem then guarantees that the distribution of subspaces orthogonal 
to k and m is (locally) integrable. I recall that (|96| ) is implied by the field equations for 
vacuum spacetimes and also for certain matter models (electromagnetic fields, ideal 
fluids, but not for Yang-Mills fields). 

In this situation, spacetime is (locally) a product manifold, M = S x T, where 
£ = R x 5*0(2) and T is perpendicular to S. Thus the metric splits 

®g = tr + g (97) 

such that <r is an invariant 2-dimensional Lorentz metric on S, depending, however, 
on y G T, and the fact that (T,g) is a 2-dimensional Riemannian space. In adapted 
coordinates 



x — t, x = (f for E; x , x for T, (9£ 



we have 



k = d t , m = d v , (99) 

and 

er = a a bdx a dx b , g = gij dx % dx\ (100) 

where a, b = 0, 1 and i, j = 2, 3. The metric functions a a b and gij depend only on the 
coordinates x % of V. 

The following functions on T have an invariant meaning 

-V:=(k\k), W:=(k\m), X := (m\m), (101) 

and we have 

cr = -Vdt 2 + 2Wdtd(p + Xd<p 2 . (102) 



We use also 



W 



p := VVX + W 2 = V^, A:=—, (103) 

in terms of which ( |102| ) takes the form 

o 2 

a = -L-dt 2 + X(dp + Adt) 2 . (104) 

It turns out that the partial trace R a a of the 4-dimensional Ricci tensor is propor- 
tional to ^Ap and p is thus a harmonic function on T, whenever R a a vanishes. This 
is of course the case for vacuum manifolds, but also for the Kerr-Newman solution, 
and in some other cases J/J . With the help of the Riemann mapping theorem one can 
then show that p is a well-defined coordinate on (T,g) (p has no critical points). It is 
then possible to introduce a second coordinate z, such that 

g = L e 2h (dp 2 + dz 2 ). (105) 

In terms of t, ip, and the Weyl coordinates p, z, assumes the Papapetrou 
parametrization 

= -fL. dt 2 + X (cbp + Adt) 2 + —{dp 2 + dz 2 ). (106) 

We emphasize once more, that the functions X, A, and h depend only on the Weyl 
coordinates p and z. 

The weak rigidity theorem tells us that on the surface on which 

i = k + Vt m) VI = -A = -W/X, (107) 

becomes null, Q is constant and £ is a Killing field on this surface, denoted by H[£] 
in what follows. In addition, H[£\ is a stationary null surface, in particular a Cauchy 
horizon. (For proofs, see M.) 



Note that f , 1-form, can be expressed as 

i = - P ^dt. (108) 

One also knows that X is well-behaved on the horizon 

X = 0(1), X~ 1 = 0{1) on#[£] (109) 

(see ||). Since (£|£) = —p 2 /X the horizon is at p = (as is well-known from the Kerr 
solution). 

Although £ is not a Killing field, the weak rigidity theorem implies that the surface 
gravity is still given by 

« a = ~(deide)L Ifl - (no) 

(A priori, the formula holds for I := k + f2# m, = 0|h[£], but dfi = on 
From (|108|) we obtain 



= - ^ dpAdt-p 2 d Q^J A tZt. 

Using the scalar products (dp\dp) = Xe~ 2h , (dt\dt) = —X/p 2 , (dp\dt) = gives thus, 
together with the zeroth law, 

k = e~ h \ H ^ = const 7^ 0. (Ill) 
The reader should verify that this gives the correct result for the Kerr solution 

T H - M 



2Mr 



H 



112) 



After these preparations, our argument proceeds similarly as in (b). The natural 
FIDO tetrad {6»} for (|D|) is (N := p 2 ) 

/N p h e h 
— dt, 6 1 = , dN, 6 2 = ^dz, 9 3 = Vx(dw + Adt). (113) 
X 2^/XN VX 



Clearly, 6° and 6 1 are again ill-defined on the horizon N = 0. We therefore pass over 
to a new basis 

P^dt+^dN, P :=d^ + n-^dN, (114) 

together with dN and dz. In order to check whether this new basis remains valid when 
the horizon is approached, we express the metric ( |106| ) in terms of it. Since 

/N - e h 
— 8 1 ^=dN, 3 = VX(9' P + A6 t ), (115) 
X 2VXN 

we find readily 

W g = -Vie 1 ) 2 + 2W 6*6* + X{6*) 2 + — 6 f dN + — dz 2 . (116) 

X X 



The determinant of the metric coefficients in this expression is and thus remains 
regular, thanks to ( fLO^j ) and ( p.lU| ). Since we postulate a regular horizon, it is then 
clear that {#*, 9 V , dN, dz} indeed forms a well-defined basis also on the horizon. 

It is now straightforward to express the Maxwell 2-form F in terms of this basis. 
Instead of ( |96| ) we obtain now 

f e h Ae h 1 

F = < Ei — — H y=(E 3 — B 2 ) > dN A 0* + E 3 VW 6 V A W 

[ 2A. 2y N J 

/AT 1 - p h 

— + B ie h A }dzA9 t + —=(B 2 - E 3 ) dNA 0* 



+B ie h dz A 6 V — {E 2 + B 3 ) dz A dN. (117) 

2^^ 



Since a := y ^ is the lapse function, we obtain once more the limiting behavior ([75|), 
and thus the basic membrane laws (]69|) and (|70|) . 



5 Magnetic Energy Extraction from a Black Hole 

As an interesting, and possibly astrophysically important application of our basic laws 
(P^) and (|oT)| ) I show now that it is possible, in principle, to extract the rotational 
energy of a BH with the help of external magnetic fields. In the next section, we will 
work out some of the details for an ideal gedanken experiment. This will serve as a 
preparation for an understanding of the Blandford-Znajek process. 

Our starting point is Faraday's induction law ( p5|) in integral form, which we write 
down once more 

EMF(C) = -|- / B+lipB. (118) 
dt J a Jc 

For stationary situations this reduces to 

EMF(C) = jipB. (119) 

In Fig. |3| we consider a stationary rotating BH in an external magnetic field (like 
in §3). The integral in (|119|) along the field lines gives no contribution and far away 
P drops rapidly (~ r -2 ). Thus, there remains only the contribution from the horizon 
(Cjj) of the path C in Fig. |: 

EMF = / i Ph B, (3 H = -n H d v = -n H u H e v (120) 

Jc H 

(only the normal component B± contributes). I recall that 0# = a(2Mrn)~ 1 , th = 
M + VM 2 -a 2 . 

We shall show in section 6 that it is possible to construct a generator such that 
(with optimal impedance matching) a maximal extraction rate equal to 



1 (EMF) 2 
4 R(C H ) 



(121) 



Figure 3: Arrangement for eq. ( |119| ). 



becomes possible, where R(Ch) is the horizon ("internal") resistance 

f di 

R{C H ) = Rh—, Rh = 377 Ohm. (122) 

Jc H 27r U 

(It should be known from electrodynamics, that we have to divide the surface resis- 
tivity Rh by the length of the cross section through which the current is flowing.) We 
shall also show that an equal amount of energy is dissipated by ohmic heating at the 
(stretched) horizon. 

Let us work this out for the special case of an axisymmetric field: Lq v B = <-> 
dig v B = 0, whence 

i g B = . (123) 

dlp 2tt V ; 

^From this we conclude that B can be expressed in terms of two potentials \l/ and g, 

£ = ^d^Ady? + (r*d(£ , (124) 

v toroidal part 

poloidal part 

both of which can be taken to be independent of (p. This is equivalent to the vector 
formulae 

B po1 = -^r V* x e v , B tor = %2 V (125) 

(Exercise). ^ is the magnetic flux function (see Fig. because the poloidal flux inside 
a tube = const} is 

" B = J = ^, *(0) = 0. (126) 

\1/ is constant along magnetic field lines, as should be clear from Fig. [|. Formally, this 
comes about as follows: 

i-dty = *(*d^ A B) = *(d^ A *B) = *(d^ A B) 
g * (d\I> A *d<f) = 0, 



(124) 




magnetic flux function 
^ = const., 



Figure 4: Axisymmetric magnetic field. The total flux inside the magnetic surface 
defines the flux function 



thus < d<H,B >= 0. 

For the closed path C in Fig. | the EMF is by ( |120D 



EMF = AV 



(N) 



n 



H 



2tt 



i.e. 



H 



2tt 



(127) 




Figure 5: Integration path for the voltage in eq. (|127| ). 



The internal resistance (|122 ) becomes 



AR H = R H 



Al 
2tt Co 



'128) 



Since, on the other hand, 



A^ = 2ti uj B± Al, 



(129) 



we find, by eliminating A/, 



AR 



H 



R 



H 



A^ 

4vr 2 Co 2 B i ' 



(130) 



Inserting ( 127 ) and (|130| ) into the expression (|121 ) for the maximal power output 
gives 



1 (AV) 2 
4 AR H 16 7T 



O 2 



(131) 



This, as well as (|127|) , have to be integrated from the pole to some point north of the 
equator (see Fig. For the exact solution in §3 we know the result for the EMF, if 
we integrate up to the equator: From (|127|) and (|67|) we get 



EMF 



2tt 



Q H An B M (r H - M) 



or (tt H = a/2Mr H ) 



EMF = a B, 



r H - M 



(r H = M+ VM 2 - a 2 ) 



(132) 




Figure 6: Plausible structure of the magnetic field close to a supermassive BH which is 
surrounded by an accretion disk (central engine of a quasar) (adapted from Ref. [1]). 

For a general situation, like in Fig. 6, we have roughly 

S A^I = ^ ~ B ± 7rr 2 H , lj 2 ~ < lj 2 > ~ — , 

and we obtain, by ( |131| ), for the power output 




(133) 



The total EMF is (V = £ AV) 



V~— Q H ^~- —B ± nr 2 H ~-(—)MB L (134) 

2tt 2tt 2Mr H H 2 VM/ V ; 

(compare this with ( |132| )). Numerically we find 

For reasonable astrophysical parameters we obtain magnetospheric voltages V ~ 10 20 
Volts and power output of the magnitude ~ 10 45 erg/s. This power is what one ob- 
serves typically in active galactic nuclei, and the voltage is comparable to the highest 
cosmic ray energies that have been detected. 

Note, however, that for a realistic astrophysical situation there is plasma outside 
the BH and it is, therefore, not clear, how the horizon voltage ( |135| ) is used in accel- 
erating particles to very high energies. 

Let us estimate at this point the characteristic magnetic field strength than can 
be expected outside a supermassive BH. A measure for this is the field strength Be 
for which the energy density B e /8tt is equal to the radiation energy density ue cor- 
responding to the Eddington luminosity 

i£= 4rfW = 13xl0 3 8 (erg/s) M 
The relation between Le and ue is 



„ 2 c 2 ( GM\ . . 

A-KT g - u E = vrr 9 cu E I r g = — ^- ] . (137) 



Thus 



J-B| = l™< (138) 
8 7r a-? r g 

giving (Mff )8 = M H /10 8 M ) 

B E = 1.2 x 10 5 M~]l 2 Gauss. (139) 

For a BH with mass ~ 10 9 M & inside an accretion disk acting as a dynamo, a char- 
acteristic field of about 1 Tesla (10 4 Gauss) is thus quite reasonable. 



6 Rotating BH as a Current Generator 

Before we come to realistic possibilities of energy extraction, we analyze in detail an 
idealized arrangement, sketched in Fig. [7]. 

We compute first the EMF around a closed path consisting of the following parts: 
We start from the equator of the (stretched) horizon along a perfectly conducting 
disk, which is supposed to rotate differentially in such a manner, that all its pieces 
are at rest relative to the FIDOs (they have thus zero angular momentum). From the 
boundary of the disk, the path continuous along a wire and through a resistive load 
(Rl) to the top of a conical conductor. Then we move down to its tip at $ = ^ 




Figure 7: An idealized current generator (adapted from [1]). 

close to the north pole of the membrane, and finally down the horizon in the poloidal 
direction to the starting point. (The conical surface is chosen instead of an infinitely 
thin wire in order to avoid divergent integrals; it replaces a wire of finite thickness.) 

The ideally conducting disk gives no contribution to the EMF. Formally this comes 
about as follows: The 4- velocity field of the disk u (whose parts move with the FIDOs) 
is e (see (|5|)) and thus —i u F = E = 0, for an ideal conductor. Using (^), this implies 
£ = 0. 

The contribution Vj, to the line integral § a£ from the load resistance far from 
the BH is what we all know from electrodynamics 

V l = IRl, (140) 

where I is the current in the wire. 

There remains the contribution Vh to the EMF along the stretched horizon 

V H = f a£ { °^ [ R H *J H . (141) 



J Ch J Ch 

The current flows along the northern hemisphere of the horizon as a surface current 
in the poloidal direction, with surface current density 

J H =^—^e d . (142) 

In order to prove that the surface current Ju has no toroidal component, we 
apply ( |119| ) to a toroidal path = const} on the stretched horizon. Along this 
ip B = -d^/2vr = 0, and thus 

a £ = =^ i 9(fi £ H = : E H ■ e v = 0. 



I— I — * 

Ohm's law in (p9|) then implies Ju • e v = 0. 
Using ([ggp in QI1I|) gives 



or 



where 



V H = IR 



HT, 



(143) 



1? 



is the total resistance of the horizon (see ( |122| )). The EMF around our closed path is 
thus 



V = V H + V L = I{R HT + R L ). (144) 
This voltage is also equal to the line integral on the right hand side in ( |119| ), 



V 



ipB. 



(145) 



This receives only contributions from the horizon and the disk (the contribution of 
the latter was overlooked in Ref. [1]). Using f3 - 
— u dip, the horizon gives 



u d v , B± = B r ¥ A = pdd, 



ipB = Q H Bj_ wh p H d-d. 
c H J$o 



(146) 



A similar contribution is obtained along the disk (Cjj) and the total voltage is given 
by 



V = £l H I B± ujh Ph d-d — / B± uj uj -j= 



I'D 



dr 



(147) 



(td is the edge of the disk). This "battery" voltageQ, and the total horizon and load 
resistances Rht and Rl determine the current I according to (|144|) . 

If the magnetic field is axisymmetric, we can use ( |123|) to write the integrand in 
( |j~45l ) as follows 



to 



148) 



The voltage is then 



V = ^*(eq) 

I 71 



edge 



UJ 



horizon 

[n H ) d V^(eq) 



(149) 



1 The part of the integral ( |145| ) from the disk is independent of the connecting path, because the 
induction law gives d(a£ — ipB) = and thus d (ipB) = inside the disk. 



where \l/(eq) denotes the value of the flux function at the equator of the horizon. The 
two pieces in the square bracket are comparable in magnitude. 

We shall see in the next section in detail how the power, dissipated as ohmic losses 



Pl = I 1 Rl (150) 

in the load, is extracted from the hole, but this is clearly at the cost of the mass of 
the BH: 

I 2 R L = —M. (151) 

Let us note that 

p ^ 2 <raF (152) 

This becomes maximal for 

Rht — Rl (impedance matching) (153) 

with 

PT aX = T~B~ ■ (154) 

This maximal extraction rate was stated in fljgjj) . 

Clearly, we can also reverse our gedanken experiment. By applying a voltage, we 
can use the BH as the rotator of an electric motor (see Fig. |8|). You see that some of 
the physics of BHs is indeed very similar to that of ordinary electric generators and 
electric motors. 



7 Conservation Laws, Increase of Entropy of a BH 

The general results which will be obtained in this section will enable us to develop a 
more profound analysis of the idealized current generator discussed above. 

Contraction of the energy- momentum tensor T^ LU with the two Killing fields d t , d v 
gives two conserved vector fields, which express the conservation laws of energy and 
angular momentum in the z-direction. We want to formulate these in the 3+1 splitted 
form. 

To this end, we consider as a preparation a 4-dimensional equation of the type 

V-J = Q (155) 

for a vector field J with source term Q. The Hodge dual of the 1-form J b has the same 
decomposition as S in (P2|): 

*J b = p + (ipp-aj) Adt. (156) 

If J = J^e^, then p is the 3-form belonging to J°, and J is the 2-form corresponding 
to j = J k 6k- Eq. (|155|) is equivalent to d (*J b ) = Qvol 4 or 

dt A [(d t — Lp)p + d(a J)\ = Q vol 4 . (157) 




Figure 8: Black hole playing the role of the rotator of an electric motor. 



Since dt = 9 u /a, the 3+1 split is 

(dt-L p )p + d(aJ)=aQvoh. (158) 
If div/3 = (as for the Kerr solution), this is equivalent to (p =: PV0I3) 

{dt-L )p + V- (a J) =aQ. (159) 

The integral form 

4 [ P = ~ [ (aj-ipp)+ [ aQvoh (160) 

generalizes the conservation law fl28"|) . We also write ( |160| ) in vector analytic form 

_ f pdV = - [ (a j - p(3) ■ dA+ [ aQdV. (161) 

Let us now apply this to the vector fields k^T^ und m fl T fJ,u (k — d t , m = d v ). 
The corresponding p and j are denoted as follows 

-k-T < — ► (e Ego , S E J, 

m-T <— (sl z JlJ. (162) 

Eq. ( [T59D gives 

(d t - Lg) e Eoo + V • (a S Eoo ) = (Poynting theorem), (163) 

(d t — Ls) Sl z + V ■ (cxSl x ) = (angular momentum conservation). (164) 



The corresponding integral formulas are 



± I £Eoo dV = - I ( a S Eoo -e Eoo f3)-dA, (165) 

d 
dt 



[ e Lz dV = -[ (aS Lz -e L j) ■ dA. (166) 

JV JdV V J 



Now we make use of @, i.e., d t = a eo — u) d v = a eo — u) Cj e v , giving us 

—k ■ T = —a e ■ T + uj m ■ T. 
Here we use the following FIDO decomposition of T: 

T = ee ®e + e ®S + S®e +T, (167) 

and obtain the relations 

£ e cc = oce + uj£l z , (energy density "at infinity"), (168) 
Seco = aS + uSL z , (energy current density "at infinity"). (169) 

Application to a Kerr BH 

The global conservation laws ( |165|) and ( |166|) are now applied to a Kerr BH. Its 
mass plays the role of the "energy at infinity" and thus, the energy conservation ( |165| ) 
with ( |169|) , gives 



dM 



/ a H S Eac -ndA = - (a 2 H S + a H Q H S L )-ndA. (170) 



dt 

Similarly, the change of the angular momentum of the BH is by ( |166|) 

-r = -\ a H S Lz -ndA- S E = d v - T • (171) 
dt J H s 

Now, we consider the entropy increase of the BH. The first law M tells us that 



_ dS H dM n dJ . , 

T °-*= (172) 

On the right hand side we insert ( |170| ) and ( |171| ), and use also (|169| ), giving us the 
generally valid formula 

dSu f 2 



T H —f L = - / a 2 H S-ndA. (173) 



Until now we have not specified the matter content. For electrodynamics we have 
S = j- E x B, and thus at the horizon a 2 H S — j- Eh x Bh (see (|78|)). In this case 
we obtain with the laws of Ampere and Ohm 

Th< ^ = --tI (E H xB H )-ndA= [ J H -E H dA= [ R H J H dA. (174) 

All of these familiarly looking expressions for the rate of entropy increase are important 
and useful. 



Let us also evaluate ( |171|) in a similar manner. First, we have 



E®E + B®B) + -[E Z + B Z ) g 



(175) 



Clearly, only the first term contributes to the integrand in ( |171| ). Using this time the 
laws of Gauss and Ampere, we find 



dJ 



a H E H + J H x B n ) ■ d v dA. 



(176) 



Note that the first term is absent if Eh has no toroidal component. 

Finally, we use (|172|) , together with (|174|) and ( |176j ), to obtain the following formula 
for the change of the mass of the BH 



dM 
~dT 



J H -E H -/3 H - [a H E H + J H x B r 



dA. 



(177) 



Application to the idealized current generator 



It is instructive to use these general results for a more detailed analysis of the 
current generator, discussed in the last section. 

We begin by computing the ohmic heating rate of the current flowing through the 
northern hemisphere (n.H.) of the BH: 



E H ■ J H dA 



n.H. 



(hi 



(142) T / 2 Z? A Q T \r (£la) Jl D 

= / / E H ^p H dd = IV H - I Rht- 



(178) 



According to the general result ( |174j ) this rate is equal to T# dSn/dt. Thus, 

T h ^ = I 2 Rht- (179) 



In order to trace the details of the energy flow, we apply the generalized Poynting 



theorem ( |165| ). For a stationary situation this reduces to 



/ (aS Eoa -e Eoo 0)-dA = O. 



(180) 



We choose the volume V such that the boundary dV consists of a horizon part Ah 
(TC S minus the inner edge of the disk), the disk Ad, and a surface Ai enclosing the 
load's resistor (see Fig. |). The second term in (|180|) does not contribute for Ah and 
Ad, and can be ignored for the load, since this is assumed to be located far from the 
horizon. Using also the relation ( |169| ) we have then 



aS E ^ ■ dA 



(a 2 S + acuS L ) -dA = 0. 



;i8i) 



dv 



± 



T 




Figure 9: The boundary dV = A H U Ad U A L for (1181 



The contribution from the horizon to the first term on the right is 



and the second term gives 



a^S* ■ n dA ^-=* Tn 



-Qh I can Sl z • nd A. 



182) 



(183) 



At first sight one might think that this is just ^H^jf (see ( |171| )). This is, however, 



not correct, because there is an additional inflow of angular momentum through the 
disk. Physically, this is clear: The external magnetic field exerts - through the Lorentz 
force on the surface current Jd of the disk - a torque on the disk, and thus on the 
BH to which the disk is locked. As in the derivation of ( |176|) we find for ( 183|) 



-Q H I a H S Lz ■ ndA 

A H 

-tt H I J H $ B± uj h 27r Cj h p H d$ 

h 



\ u H (J H x B n ) ■ e 9 dA 



(142) 



-IVLh I B±u H pHd'8, 

'1?0 



where use has been made of the fact, that Eh has no toroidal component. Together 
with ( |146| ), we obtain for this horizon contribution 



-VL H I a H S Lz ■ ndA 

I Ah 



-I I ipB. 

>Ch 



;i84) 



Note that the integral on the right is not the total voltage (|145|) . However, we shall 
see shortly that the disk contributes the remaining part of —IV. 

Only the last term in (|181|) gets contributions from the disk (since E = in the 
disk). Using the expression ( |175| ) for Sl b , this becomes 



a to Sl z ■ n dA - 

Ad J Cd 



a uj — du, ■ B 2ttu dl. 

47T V 



;i85) 



— * — * 

Between the parallel component Bn and the surface current Jd of the disk we have 
Ampere's relation of ordinary electrodynamics: 



aB\\ = Ait J D x n = An - — - e r x n. (186) 



Making use of this, ( |185| ) becomes 

- J auoS Lz -ndA = -I J (3 ■ (e r x B^j dl = I J $ ■ (j} x B^j ■ dl 



or 



auS Lz -ndA = - i p B, (187) 

A D Jc d 



as already announced. 

Finally, the contribution to (|181|) is what we are used to: 



/ aS Eoo -dl= [ -*-(£ x B)dA = I 2 R L . 
Ja l ~ J a l 47r 

All together, the Poynting theorem (|181|) gives 



horizon 




I Rth-I j ipB-I j ipB = —I Rl- (189) 



load 



^From the derivation it is clear how to interpret this result. First, we note that 
the left hand side is the energy (measured at infinity) which flows per unit time into 
the BH, and thus is equal to dM/dt. According to ( |179| ) the first term is T# dSn/ 'dt 
(ohmic dissipation). The second term in ( |189|) is that part of Qh dJ /dt which is due to 
the torque acting on the surface current density Ju (see ( |171| ) and ( |176| )). The third 



term resulted from the disk and gives an additional contribution to the change of the 
rotational energy, which flows through the disk into the BH. For clarification we note 
that for a closed surface A surrounding the disk we have 

a Se^ ■ dA = = (b auj Sl z • dA, 
A Ja 

since S vanishes in the disk. Therefore, the rotational energy which flows through the 
inner edge of the disk into the BH is (see ( |171|) ) 

Qh / ot- Sl z ■ dA — — I auo Sl z ■ ndA ^ — ^ — I ipB. 
J+-\ Ja d Jc d 

The total change of the rotational energy of the BH thus is 

^h-^ = -IfipB = -IV 



c 



(144) 



-I\R HT + R L ). (190) 



The presence of the disk made the discussion a bit complicated, but it is nice to 
see how the various pieces combine. Schematically, we have 

dM _ rp dSH I n dJ 

\\ || || ;> (191) 

—I Rl I Rht —I (Rht + Rl)- 

The energy (at oo) dM/dt flows - partly through the disk - down the hole and is, of 
course, negative because the spin down overcomes the ohmic dissipation. This part of 
the energy balance results, therefore, in an outflowing energy which is carried without 
loss by the electromagnetic fields to the load resistors, where it is dissipated at the 
rate I 2 Rl- Netto, we obtain, of course, 

8 Blandford-Znajek Process 

Let us return to Fig. [| in which a plausible magnetic field structure around a super- 
massive BH in the center of an active galaxy is sketched. Rotation and turbulence in 
an accretion disk can generate magnetic fields of the order 1 tesla, as we estimated at 
the end of section 5. 

Close to the BH one expects a force-free electron-positron plasma, which comes 
about as follows. Imagine first that there are no charged particles in the neighborhood 
of the hole. In this case unipolar induction generates a quadrupole-like electric field 
similar to that we found in section 3 (see Fig. |2]). Close to the BH the magnitude 
of this electric field is E ~ Ba/M ~ 3 x 10 6 (Volt/cm) (a/M). Between the horizon 
and a few gravitational radii along the magnetic field lines this gives rise to a voltage 
V ~ Er H ~ B a ~ 10 20 Volt for a BH with mass M ~ 1O 9 M (see section 5). 
In this enormous potential stray electrons from the disk or interstellar space will be 
accelerated along magnetic field lines to ultrarelativistic energies. Inverse Compton 
scattering with soft photons from the accretion disk leads to 7-quanta which in turn 
annihilate with soft photons from the disk into electron-positron pairs. These will again 
be accelerated and by repetition an electron-positron plasma is generated which can 
become dense enough to annihilate the component of E along B. The electric field 
is then nearly orthogonal to B, up to a sufficiently large component which produces 
occasional electron-positron sparks in order to fill the magnetosphere with plasma. 
(Similar processes are important in the magnetospheres of pulsars. It is likely that all 
active pulsars have electron-positron winds.) 

Fields with E ■ B = are called degenerate. We assume in what follows that in 
the neighborhood of the BH, where the S-field is strong, the electromagnetic field is 
force-free, which means that the ideal MHD condition 

p e E + jxB = (193) 

is satisfied. Clearly, in a force-free plasma E is perpendicular to B. Furthermore, E has 
no toroidal component for an axisymmetric stationary situation. This is an immediate 
consequence of the induction law: Applying ( 25) for a stationary and axisymmetric 



configuration to the path C in Fig [10], we obtain 



C JC 2 7T j c 



Similarly, Ampere's law in integral form (|2^j) reduces to 




= const. 



Figure 10: Path for O and fll95p . 



aH = 4:TT <b aj = 4irl, (195) 

C J A 

where / is the total upward current through a surface A bounded by C. This gives for 
the toroidal component of the magnetic field 

^ tor= 2J 



Using this in ( |1 24| ) gives 



1 21 
B = — dvfAd^9 + — *d(p. (197) 
2 7r / v a 

poloidal toroidal 

For later use, we also note the following: The continuity equation ( p4|) reduces to 

d(aj) = 0. (198) 
Moreover, Lg (a J) = 0, we have d (ig a J) = 0, thus 

ia v (<xJ) = 7rU (199) 

or 

ajp 01 = dJ Adyp. (200) 

2 7T 

Clearly, the potential I is the current in ( |195| ). 

Because S is poloidal, we can represent the electric field as follows 

£ = i$ F B (E — -v F x B), (201) 

where Vp is toroidal. Let us set 

vp —: — (Tip — uj) uj e v . (202) 
a 



For the interpretation of flp note the following: For an observer, rotating with angular 
velocity Q, the 4-velocity is u = it* (pt + fi d^). On the other hand, u = 7 (eo+v), where 
v is the 3- velocity relative to a FIDO. Using also d t = a eo + (3 we get Vtd v = av — (3 
or 

u = — (Q — a;) d v = — (f2 — u) £0 e v . (203) 
a a 



This has the same form as (P02j) and, therefore, ilj? is the angular velocity of the 
magnetic field lines. 

Next, we show that Qp is only a function of The induction law gives 

d(a£) = LpB = —L u , dv B = —uLg lp B — du}Aia v B 

= — duAdty. 

On the other hand, ( M) and fl20p give 

d(e*£) = d(ai^B) = d((n J ,-o;)i^B) = -dr^^d*J 

— A dtf . 

2tt / 

By comparison, we get dQp A d^ = ==>- Qf = Qp(^>). The calculation above also 
shows 

aS = -^^dV, (204) 

2 7T 

i.e., is perpendicular to the surfaces = const}. 
Specializing to the horizon gives 

E H = -(tt F - tt H ) Qe v xB ± . (205) 



The representations ( p.97| ) of B, and ( |2U4j ) for £ will be important in the final section 
9. 

Now, we proceed as earlier in section 5 and consider again the closed path C in 

Fig. 

We already had the relations 



AV = — n H AV, (206) 

2 7T 



and 



AR « = R h a 2 ~2 p 207 ) 

4:7T Z U! Z B± 

(see ( 127|) and ( |130|) ). As in ([143 ) we obtain for the horizon voltage 

AV H = IAR H . (208) 
This contribution can alternatively be computed with ([204 ): 

AV H = [ aS= Qh ~ Qf Av&, (209) 



Figure 11: Magnetosphere of a rotating BH (adapted from [1]). 



The total voltage AV of C is the sum of AVh and the voltage drop AVl in 
astrophysical load region (see Fig. |TT| ). For the latter we obtain from (|204|) ( 
io ~ 0) 

AV L = ^n F AV 

and, of course, also 

AV L = I AR L , 

where ARl is the resistance of the load region. From this equations, and 

= AV H + AV L , 

we immediately find the relations 

AV L V f AR L 

av h ~ n H -n F ~ AR H ' 

i = - — — = -(n H -n F )ti 2 B ± . 

AR h + ARl 2 v t! 
The ohmic dissipation at the horizon is as in ( |1 79|) 

T H d -^ = I*AR H = IAV H 
at 

4tt 

The power APl deposited in the load is thus 

AP L = J 2 ARl = I AV L ^ ^ - Qf) ^ B ± A*. 



As in (|153|) , AP L becomes maximal for AR H = AR L , which is the case for (see ( |213| )) 



n F = ^n H . (217) 

Then the estimates ( |133| ) and ( |135| ) hold. 

Whether the crucial condition (|217|) is approximately satisfied in realistic astro- 
physical scenarios is a difficult problem for model builders. Therefore, the question 
remains open whether the process proposed by Blandford and Znajek || is important 
for the energy production in active galactic nuclei. It could, however, play a significant 
role in the formation of relativistic jets. 

In this connection it should be mentioned that the Blandford- Znajek process may 
play an important role in gamma-ray bursts, as has been suggested, for instance by 
Meszaros and Rees p0 |. One of the proposed models involves the toroidal debris from 



a disrupted neutron star orbiting around a BH. This debris may contain a strong 
magnetic field, perhaps amplified by differential rotation, and an axial magnetically- 
dominated wind may be generated along the rotation axis which would contain little 
baryon contamination. Energized by the BH via the Blandford- Znajek process, a nar- 
row channel Poynting-dominated outflow with little baryon loading could be formed. 
The latter property is crucial for explaining the efficient radiation as gamma rays. 
(For a review of this topic, see for instance [|TT|].) 

9 Axisymmetric, Stationary Electrodynamics of 
Force-Free Fields 

In this concluding section we discuss the electrodynamics of force-free fields around 
a BH a bit more systematically. The main goal is to derive the general relativistic 
Grad-Shafranov equation, whose nonrelativistic limit plays an important role in the 
electrodynamics of pulsars. 

For stationary fields, Maxwell's equations (^) become 

dB = 0, d(a£)=LpB, (218) 
d*£ = 4vrp, d(a*B) = Arc a J -Lp*£. (219) 

^From now on we assume also axisymmetry. Then Lp p = d ip p = 0, and hence the 
continuity equation ( |24"D reduces to 



d(aj) = 0. (220) 



I recall the representation ( |197|) of B in terms of the potential \1/ and /: 



1 21 
B= — d^Adcp + — *dip. (221) 
2 7r a 

and / are independent of (p; their physical significance has already been discussed: 
Using the notation in Fig. [K] we have 

^ = [ B, 1= [ a J. (222) 

J A J A 



At this point we make the simplifying assumption that the ideal MHD condition 
( |193j ) holds everywhere. The fields are thus force-free, and we know from section 8 
that the £-field is poloidal and can be represented as (see ( |204|) ) 



a £ 



Vtp — 00 
2tt 



The representation ( |200|) of the poloidal current 

1 



a J 



pol 



2tt 



dl A d^. 



(223) 



(224) 



was obtained without assuming the ideal MFD condition. But if this is assumed, the 
toroidal part of j x B in ( |193|) has to vanish, which means that j7" po1 and £> po1 are 
proportional to each other. Comparison of Q221| ) and (|224|) shows that dJ must then 
be proportional to d\l/, thus / is a function of \P alone: 



/ = /(*). 

We have shown already in section 8 that fl F is also a function of 

n F = n F (v). 

According to ( p21| ) and ( 224) we now have 

1 dl 



jpol 



a d$> 



^gpol 



(225) 



(226) 



(227) 



The total current j can now be represented as 

1 dl s 
j = p e v F + - — B. 

a dw 



(228) 



In view of ( |227|) , the poloidal part of this equation is certainly correct, and the toroidal 
part on the right hand side is chosen such that j x B = p e v F x B = — p e E (see (|201|) ), 
which is just the ideal MFD condition. 

Our goal is now to derive - for given functions (|225|) and ( 226 ) - a partial differ- 
ential equation for the potential We shall achieve this by computing the toroidal 
part of the current in two independent ways. First, we take the toroidal part of (|228|) 
and obtain with B tOT = M- (see 



J 



tor 



Pe 



>-lu 1 dl 21 

u + - — — . 

a a aW aoo 



Here, we also eliminate p e : From (|223|) and ( p!9| ) we deduce 



1tt 2 p 



-d * 



d^ 



a 



(229) 



i.e., 



ln 2 p e 



-V 



fip — OJ 



a 



(230) 



We thus arrive at 



J 



tor 



1 Qp — UJ 



■ujV 



Vtp — uj 



a 



21 dl 



a 2 uj dW 



(231) 



On the other hand, this quantity can also be obtained from Maxwell's equation 
191): 



(232) 



A7iaJ tor = [d(a * B)] tor + [L p * S] tor . 
For the first term on the right we have 

2tt [d(a * B)] tOT = 2 vr d(a * B tOT ) © d [a * (dtf A d^)]. 



(233) 



Now, we use the following useful general identity, whose proof is left as an exercise: 
Let x be a p-form and m a Killing field. Assume also x = 0, then 



5 (m A x) — — m A 8 X? 

where 5 is the co differential. 

For x — i a /& 2 ) d^, m = 3^, m = uo 2 dip this gives 



(234) 



or 



Thus, (12331) becomes 



5(ad$A dy?) = 5 (J^ d^ to 2 dip 
d*[a d\I> A dip] = S (J^ d^ Co 2 * dip. 



27rd(a * B p ' 



oh 



a 



-V • ( — WJ uj * dip. 



(235) 



Next, we turn to the second term in ( 232|) . By (|223|) we have 

Qp — uj 



-La * 8 = Lr 



2rra 



* d\I> 



If x denotes the square bracket, we can write Lp x = L^q x — —UJ Lq v x~ du Aig X-, 
and obtain 



*Lfj * S 



Hp — UJ 



»v w (d* AO' dip) 



2na v 

With (|35D and (|36|) we obtain for the Hodge-dual of (|232|) 

f2p — UJ 



(236) 



4ira * J 



tor 



1 
2~K 



— V- (^W) Co 2 dip 



UJ Z 



2ira 



Vuj ■ uj dip. 



Since *J tOT = j tOT 4 cD 2 dy>, we get 



47ra; 
— —j 

UJ 



tor 



-— v- — w 



2tt 



a 



cu 2 / 2%a 



(237) 



Inserting Vu = V(u> — Q F ) + (dQ F /dty) V\I/ leads finally to our second formula for 

„' tor . 



87r 2 j tor = --v-(^v*) + ^(n F -u)v*-v(n F -u) 

-^(^-^)^(W) 2 . (238) 
cr dw 

Comparison of ( |232| ) with ( |238[ ) gives, after a few steps, the following generalized 
Grad-Shafranov equation: 



V 



a 
ti 2 



a 2 



, fir — uj dVLp , -* n9 16 7r 2 (if 

+ — -f- W 2 + ^— J — = 0. 239 

' a d^> y ' auu 2 d^> ' 



The integration of the original equations (for axisymmetric stationary situations) is 
reduced to this single partial differential equation. and f2p(\I/) are free functions^, 
and if \I/ is a solution of ( |239| ) the electromagnetic fields are given by (|221|) and ( |223| ), 



while the charge and current distribution can be obtained from Q230|) , (|224j) , and 



A limiting case of (|239|) is known from the electrodynamics of pulsars: Ignoring 
the curvature of spacetime and using that Q F is equal to the angular velocity Q of 
the neutron star (to derived shortly), we get the pulsar equation: 



-fl r „ 91 . 1 16 7T 2 dl 
V ■ <^ — 1 - Vi 2 Cu 2 W \ + —— I — = 
[a; 2 L J J on 2 



(240) 



(a) is the radial cylindrical coordinate in flat space). 

The equation Q239| ) holds also outside an aligned pulsar, since all the basic equa- 
tions in section 2 remain valid there. However, the different boundary condition at 
the surface of the neutron star implies Q F = Q, as we now show. In the interior of the 
neutron star the 3- velocity is (see ( |203| ) ) 

v = -(n-u)B v . (241) 
a 

Since the neutron star matter is ideally conducting, the electric field there is 

£ = ijB = —(tt-uj)dty, 

2ira 



if the £>-field is poloidal (first term in ( 221|) ). At the boundary this has to agree with 
( |223j ), implying Cl F = Cl. 

Another remark should be made at this point about the interior of the neutron 
star. We showed in section 8 that dQ F A = implying df2 A d\I/ = inside the 
star. It will, in general, not be possible to represent $ as a function of Q, since \& has 
to satisfy the Grad-Shafranov equation for = inside the star. Therefore, the 
rotation must be rigid, Q = const. 



2 These are, of course, restricted by boundary conditions, but we do not discuss this here. 
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